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FOREWORD 


This  report  is  an  advance  copy  of  a paper  that  will  appear  in  the 
Journal  of  Engineering  Mathematics  in  late  1977.  The  report  is  being 
issued  so  the  information  can  be  disseminated  before  the  journal 
publication  date. 

ABSTRACT 

Expressions  are  derived  for  the  two-dimensional  surface  elevation 
resulting  from  an  oscillatory  translating  surface  pressure  distribution. 
The  surface  elevation  is  given  a-»  the  sum  of  tour  terms,  each  of  which 
is  associated  with  an  improper  integral  having  a simple  pole  singularity. 
Results  are  presented  for  the  delta  function  and  the  uniform  spatial 
pressure  distribution. 

The  mean  work  done  on  the  fluid  per  unit  time  by  the  delta  function 
pressure  distribution  is  given.  Numerical  results  are  presented  for  the 
surface  elevation  resulting  from  the  uniform  pressure  distribution. 

ADMINISTRATIVE  INFORMATION 

This  work  was  funded  by  the  Amphibious  Assault  Landing  Craft  Program 
of  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center, 

Task  Area  SI 41 7 , Task  14174,  Work  Unit  Number  1-1180-004. 
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INTRODUCTION 


The  current  development  of  air  cushion  supported  marine  vehicles 
has  kindled  a renewed  interest  in  the  analysis  of  water  waves  produced 
by  surface  pressure  distributions.  Initially  the  primary  application 
was  for  the  prediction  of  wavemaking  drag  produced  by  translating  pres- 
sure distributions  having  various  planform  shapes.  More  recently, 
unsteady  problems  have  been  investigated  in  order  to  obtain  insight 
into  various  aspects  of  the  dynamic  performance  of  air  cushion  supported 
vehicles  [1],  [2].* 

Stoker,  [3]  Kaplan,  [4]  and  Wu  [5]  treated  the  two-dimensional 
problem  for  a harmonic,  uniformly  translating  delta  function  pressure 
distribution.  Kaplan  obtained  expressions  for  the  surface  elevation 
in  the  far  field  while  Wu  obtained  asymptotic  results  for  the  velocity 
potential  and  surface  elevation  in  both  the  near  and  far  field.  Stoker 
discussed  the  qualitative  behavior  of  the  solution,  and  obtained  results 
for  the  near  and  far  field  for  the  zero  speed  case. 

Debnath  and  Rosenblatt  [6]  treated  the  two-dimensional  finite  depth 
problem  using  generalized  function  theory  to  obtain  an  asymptotic  solution 
The  same  technique  has  been  applied  recently  by  Pramanik  [7]  to  the 
two-layer  fluid. 

Lighthill  [8]  analyzed  the  qualitative  nature  of  the  three-dimensional 
wave  pattern  produced  by  an  unsteady  translating  pressure  distribution. 
More  recently,  Tayler  and  Van  Den  Driessche  [9]  used  ray  theory  to  obtain 
qualitative  results  for  the  three-dimensional  wave  pattern  produced  by  a 

♦References  are  listed  on  page  21. 
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periodic  translating  submerged  source.  The  problem  treated  in  this  paper 
is  closely  related  to  the  translating  submerged  source  of  pulsating  strength. 

A useful  survey  and  discussion  of  the  literature  on  the  pulsating  source 
has  been  compiled  by  Wehausen  and  Lai  tone  [10]. 

This  paper  presents,  for  the  first  time,  an  entirely  analytical  result 
for  the  two-dimensional  surface  elevation  that  is  valid  for  the  entire  field. 
Results  are  given  for  both  a delta  function  and  a uniform  pressure  distribution. 
It  was  possible  to  obtain  an  analytical  solution  to  the  problem  by  directly 
evaluating  improper  integrals  arising  from  a Fourier  spatial  transform.  The 
integrals  were  evaluated  by  a careful  examination  of  the  singularities,  algebraic 
manipulation  and  proper  choice  of  contour  of  integration.  The  results  are  given 
in  terms  of  well-known  transcendental  functions.  The  results  will  be  extended 
to  three  dimensions  in  the  future,  but  a detailed  analysis  of  the  two-dimensional 
dispersion  was  considered  necessary  to  clarify  the  fundamental  nature  of  the 
wave  field.  The  methodology  used  is  equally  applicable  to  the  three-dimensional 
problem. 

The  two-dimensional  problem  is  treated  here  for  an  i rrotational , 
incompressible  invlscid  fluid  of  infinite  depth  using  linearized  potential 
theory.  The  formulation  of  the  problem  generally  follows  Doctors'  three 
dimensional  analysis  [2]  and  also  parallels  Wu's  [5]  two-dimensional  analysis. 

The  major  departure  from  Doctors  and  Wu  comes  in  the  evaluation  of  the 
integral  forms  by  contour  integration  and  the  application  of  Cauchy's  theorem. 


3 


DERIVATION  OF  VELOCITY  POTENTIAL  AND  SURFACE  ELEVATION 


The  coordinate  sysem  and  notation  are  shown  in  Figure  1.  The  surface 
pressure  distribution  p(x,t)  is  translating  to  the  right  with  a speed  c 
relative  to  the  fixed  coordinate  system.  The  surface  elevation  is  denoted 
as  z(x,t).  The  fluid  is  considered  to  be  irrotational  and  incompressible 
so  that  a velocity  potential  <j>(x,z,t)  exists.  The  potential  obeys  the 
following  relation: 

a2  <t>(x,z,t)  = 0,  (1) 

where  is  the  two-dimensional  Laplacian  operator.  The  fluid  velocity 
is  the  gradient  of  the  potential,  or 

u = <i> 

* (2) 

W = 

where  u is  the  x component  and  w is  the  z component  of  the  velocity, 
and  the  subscripts  x and  z represent  partial  di fferentation  in  each 
respective  direction. 

The  kinematic  free  surface  boundary  condition  is  written  in  linearized 
form  as  follows : 

(O  + cZ  - Z.  = 0 (3) 

z= 0 xi 

The  linearized  dynamic  boundary  condition  is: 

(<t>t  - c<t>x  + u<j>)  z=0  = •(£  + gZ),  (4) 
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where  p is  the  Rayleigh  viscosity.  The  temporary  introduction  of  Rayleigh 


viscosity  proves  to  be  useful  in  the  interpretation  of  improper  integrals 

★ 

derived  later  in  the  paper.  The  combined  free  surface  boundary  condition 
may  be  written  as 

[*tt  • 2c*xt  + + g*Z  + u($t  - c*x)]z=0  = 'p  (pt  • cpx}-  (5) 

One  may  also  write  the  following  condition  for  the  fluid  of  infinite  depth: 

<t>z  = 0,  z ■+  - 00  . (&) 

Since  the  case  of  harmonic  time  dependence  is  being  treated,  one  may  write 
for  the  pressure 

p(x,t)  = p (x)e"1ot,  (7) 

where  a is  the  radian  frequency. 

It  is  convenient  to  solve  for  the  response  due  to  a delta  function 
pressure  distribution  p(x)  = 6(x).  The  resulting  surface  elevation  will  be 
denoted  c(x,t)  and  the  velocity  potential  $(x,z, t).  The  response  to  an 
arbitrary  spatial  pressure  distribution  p(x)  may  be  obtained  from  c and  $ 
using  superposition  integrals  as  follows: 


* The  fictitious  Rayleigh  viscosity  is  frequently  used  as  a mathematical 
artifice  to  shift  pole  singularities  off  the  real  axis.  This  enables  one 
to  properly  interpret  the  integration  path  after  the  fictitious  viscosity 
is  removed.  (See  Wehausen  and  Laitone  [10],  p.  479).  The  problem  could 
have  been  formulated  using  the  Navier-Stokes  equations  and  modified  boundary 
conditions  similar  to  Miles'  treatment  of  the  Cauchy-Poisson  problem  [11]. 

A treatment  of  this  sort  would  be  more  satisfactory  from  a physical  stand- 
point. Either  way,  one  obtains  the  same  result  after  suppressing  the  viscosity 
The  introduction  of  a weak  internal  damping  mechanism  is  not  uncomnon  in 
other  areas  of  classical  physics.  For  example,  the  author  in  Reference  [12] 
has  applied  Voigt  viscoelasticity  to  the  problem  of  acoustic  reflection 
from  a solid  halfspace.  The  dissipative  mechanism  facilitated  the  inter- 
pretation of  a spatial  Fourier  integral  form. 
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0 ( x , z , t ) = / p(0  0(x-c#z,t)dc 


and 


Z(x 


L 

,t)  = f p(e)  t(x-?;,t)d5 

— oo 


(8) 


One  may  formally  express  the  velocity  potential  and  surface  elevation 
in  a Fourier  integral  form: 

*(x,z,t)  = f (k,z,t)eikxdk  (9) 

J — CD 


and 


c(x 


f oo 

,t)  = / c’(k,t)elkxdk, 

> _ oo 


(10) 


where  0'  and  c'  are  the  transformed  potential  and  surface  elevation. 

From  (1)  and  (6),  one  sees  that  the  transformed  potential  may  be  written 


0'  = A(k,t)e 


k|? 


ni; 


As  one  is  concerned  only  with  the  steady  state  response,  the  time 
dependence  of  the  velocity  potential  and  surface  elevation  will  be  of 
the  same  form  as  the  applied  pressure  given  in  equation  (7).  In  this 
case  the  time  differentiation  indicated  in  the  boundary  conditions  reduces 
to  multiplication  by  -io. 

The  transformed  surface  elevation  may  then  be  expressed  as  follows 
from  equations  (3)  and  (11): 

c‘(k,t)  = (12) 

One  solves  for  the  potential  function  A(k,t)  by  substituting  (11)  and 

(7)  into  the  transformed  form  of  the  combined  free  surface  boundary  condition 

(5),  giving  as  a result 


A(k , t) 


■iot 


i (kc  + a) 

{ |k|g  - u2/4  - (7+  kc  + i u/2)2  1 


(13) 
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From  (12)  and  (13),  the  transformed  surface  elevation  is: 


c‘(k,t)  = 


{ | k | g - m /4  (a  + kc  + i u/2)  } 


To  eliminate  the  absolute  signs  in  equation  (14)  one  may  break  the 
integral  form  (10)  into  two  regions,  one  along  the  negative  real  axis 
and  the  other  along  the  positive  real  axis.  After  some  manipulation  the 
following  result  is  obtained: 


c(x,t)  = 


elkxkdk 


e’lkxkdk 


J2Ik,o,c 


D-,  = ( a + kcr  + i p(o  + kc)  - kg 


D o = ( a - kcr  + i u(a  - kc)  - kg. 


EVALUATION  OF  INTEGRALS  BY  CONTOUR  INTEGRATION 
To  evaluate  the  integrals  in  equation  (15)  one  must  first  determine  the 
nature  of  the  singularities  of  the  integrands.  No  branch  point  singularities 
are  evident,  but  the  two  denominators  D-j  and  each  have  two  distinct  zeros 
corresponding  to  simple  pole  singularities.  The  second  denominator  may  be 
written  as  follows: 


Do (k , a,c)  = cZ(k  - Kj  (k  - K9), 
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where 


= Kq(2  + a - 2 /I  + a) , 
l<2=Ko(2  + a + 2 /!  + a) , 

K0  = g/ (4c2 ) » 

a = ao^  + 1 e)  * 

_ 4c o 
ao  " g 

e = y/(2o). 


and 


The  damping  is  taken  to  be  small,  so  e « 1. 

The  migrations  of  K-j  and  as  the  frequency  parameter  aQ  increases  are 
shown  in  Figure  2a.  One  sees  that  both  poles  remain  near  the  real  axis 

for  all  values  of  a . 


One  may  write  for  the  first  denominator: 


where 


and 


D«|  (k,a,c)  = c ( k - K3)  (k  - K4), 
K3  = Kq  ( 2 - a - 2 /T^-), 

K4  = K0  ( 2 - a + 2 /T^T) . 


(17) 


The  migration  of  these  poles  as  a increases  is  shown  in  Figure  2b.  The 
poles  remain  near  the  real  axis  until  a approaches  unity,  at  which  point 
both  poles  begin  to  move  away  from  the  real  axis. 

The  expression  for  the  wave  elevation  may  be  simplified  by  substituting 
(16)  and  (17)  into  (15)  and  applying  a partial  fraction  expansion  to  each 
term.  The  resulting  expression  for  the  surface  elevation  is 
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dx,t)  = 


Kl 

1 Ux, 

/VU  1 


K.)  + — Io(x,Kj 

1 * * 


K K 

-3—  I3(x.K3)  + — I4(x,K4)  , 
1- a 1" a 


where 


I1(x,K1)  = 


Jo  K - K1 


Iz(x,K2)  = 


K - K- 


I x , 1(3) 


K - K. 


^ 4 ( x » ^4 ) = 


elk*  dk. 
K - K, 


The  problem  now  reduces  to  evaluation  of  the  improper  integrals  I-j , 
I2»  I35  and  given  in  (18).  At  this  point  the  artificial  internal 
dissipation  may  be  eliminated  by  setting  the  Stokes  viscosity  coefficient 
\i  to  zero.  The  dissipation  was  introduced  to  determine  in  which  quadrant 
the  poles  , K^,  K^,  and  lie.  When  the  Stokes  viscosity  u is  set  to 
zero  the  nondimensional  frequency  parameter  a in  equation  (16)  becomes 


4oc 

a = a = — — 

o g 


The  poles  K3  and  now  lie  on  the  real  axis  for  a < 1 and  Kj  and  K2  are 
real  for  all  a.  For  a > 1,  the  complex  poles  K3  and  may  be  expressed 
in  exponential  form  as  follows: 
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aK  e 
0 


K 


3 


where 


K4  = aK0e_1>  , (a  > 1), 


\\)  = tan 


2v/o^T 

2-a 


(20) 


The  paths  of  integration  for  the  improper  integrals  in  (18)  must  be 
reinterpreted  when  the  artifical  dissipation  is  eliminated.  The  K-j  and 
poles  for  all  values  of  a and  the  and  poles  for  a < 1 now  lie 
on  the  real  axis,  so  the  paths  must  be  indented.  Each  path  is  idented  so 
that  the  pole  lies  on  the  same  side  of  the  integration  path  as  it  did  when 
dissipation  was  present.  The  expressions  for  the  integrals  may  be  written 
as  follows  for  the  nondissipati ve  medium: 


f -ikx 

Il(x’Kl)  = r J Tmc,  dk 


I2(x,k2) 


I3(x,K3) 


I4(x>K4)  = 


(21) 


where  the  paths  r-j,  r^,  and  are  now  in  the  complex  k-plane  as 
indicated  in  Figure  3. 
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SECURITY  CLASSIFICATION  OF  THIS  PACE  (Wh*n  D.r.  Bnffd) 


The  following  symmetries  exist  between  the  integrals  in  equation 
(21).  These  follow  from  the  paths  shown  in  Figure  3 and  equation  (21). 


I2(x,K)  = I^xJC), 
I3(x,K)  = I^-x.K) 

I4(x,K)  = I*(-x,K), 


(22) 


where  the  asterisk  denotes  a complex  conjugate  and  K is  real.  The  symmetry 
relations  reduce  the  number  of  integrations  from  eight  to  two. 

One  starts  by  integrating  I-j,  as  given  in  equation  (21).  The  behavior 
of  the  exponential  term  is  exploited.  One  notes  that  for  positive  x the 
exponential  of  the  integrand  for  Ij  has  a negative  real  part  in  the  lower 
half  of  the  complex  k plane.  For  x > 0,  1-j  is  evaluated  by  closing  a contour 


in  the  fourth  quadrant  as  shown  in  Figure  4a.  The  pole  is  excluded  from 
the  contour  because  of  the  indentation.  After  applying  Cauchy's  integral 
theorem  around  the  closed  contour,  one  has 


I 


+ 

1 


(23) 


where  the  superscript  plus  sign  denotes  the  solution  on  the  positive  x 
axis.  After  some  manipulation,  the  integral  I^+  may  be  expressed  in  terms 
of  auxiliary  exponential  integrals  [13]  as  follows: 

l/  (x,  Kj)  = g(KjX)  + i* f ( K -j x ) . (24) 
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One  evaluates  1^  for  x < 0 (denoted  1^")  in  a similar  fashion.  The 
contour  is  closed  in  the  first  quadrant  as  shown  in  Figure  4b.  Noting  that 
the  pole  is  now  inside  the  contour,  one  may  apply  Cauchy's  residue  theorem, 
giving 


!i  + 


-ikx  .. 
e dk 

k-K, 


2 ie-iKlx. 


(25) 


After  similar  manipulation,  the  expression  for  1^’  reduces  to 


Ij~(x,Kj)  = g(-K^x)  + i[2TTe  "Klx  - f ( -K-j x ) ] , x > 0.  (26) 


Because  of  symmetry  (22),  the  expressions  for  I^+  and  may  be  obtained 
from  (24)  and  (26)  by  substituting  for  K-j . 

The  second  symmetry  relation  (22)  is  exploited  to  obtain  the  following 
expressions  for  and  I^~  for  a < 1: 

I3+  (z.K3)  = g(K3x)  - i[f(K3x)  - 2*e1K3x]  (27) 


and 

I3"(x,  K3)  = g(-K3x)  + if(-K3x). 


(28) 


Finally,  the  third  symmetry  relation  is  used  to  evaluate  and 
1 4 for  a < 1 : 

I4+(x,K4)  = g(K4x)  - i f(K4x)  (29) 

and 

I4‘(x,K4)  = g(-K4x)  + 1[f(-K4x)  - 2*eiK4x].  (30) 
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To  evaluate  and  1^  for  a > 1 the  locations  of  the  complex  poles 
and  must  be  taken  into  account.  From  Figures  2b  and  3e  and  equation 
(20)  one  sees  that  lies  in  the  upper  half  plane  and  in  the  lower 
half  plane.  Both  poles  lie  in  the  right  half  plane  for  a < 2 and  in  the 
left  half  plane  for  a > 2. 

For  a < 2 the  expressions  for  1^  and  I4  given  in  (27),  (28),  (29)  and 
(30)  apply  because  the  poles  still  lie  in  the  same  quadrants  as  they  did 
*or  a < 1.  However,  the  arguments  of  the  auxiliary  exponential  integrals 
become  complex  and  the  behavior  of  the  residue  terms  changes  because  the 
poles  and  are  complex.  The  residue  terms  in  equations  (27)  and  (30) 
for  I and  1^",  instead  of  representing  unattenuated  surface  waves,  now 
have  exponential  attenation  as  one  moves  away  from  the  disturbance. 

For  a > 2 the  and  poles  both  lie  in  the  left  half  plane.  No 
residue  terms  occur  in  this  case,  as  neither  pole  lies  inside  the  integration 
contour.  Therefore,  both  1^  and  1^  consist  solely  of  auxiliary  exponential 
integral  terms  with  complex  arguments.  Equations  (27)  and  (30)  still  apply 
if  the  residue  terms  are  dropped  in  the  expressions  for  I^+  and  1^  . 

MEAN  WORK  RATE  AND  RADIATION  OF  ENERGY  INTO  THE  FAR  FIELD 
The  rate  at  which  energy  is  carried  away  from  the  pressure  disturbance 
by  the  free  wave  system  is  of  importance  because  of  its  association  with  the 
work  done  by  the  pressure  on  the  fluid  in  the  near  field.  In  the  far  field 
the  disturbance  consists  of  four  surface  waves  for  0 < a < 1 and  two  for 
a > 1.  The  energy  efflux  may  be  written  following  Lamb  [14],  as: 

4 

^OUT  = Z cg  En'  (31) 

n=l  n 


L 
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where  is  the  mean  energy  efflux  through  the  boundary  of  a control  volume 
moving  with  the  disturbance,  c is  the  group  velocity  of  the  nth  wave  in  the 
moving  coordinate  system  and  En  is  the  mean  energy  per  unit  surface  area.  In 
equation  (31)  the  relative  group  velocity  is  taken  as  positive  when  directed 
away  from  the  disturbance.  Calculating  the  group  velocities  from  the  wave- 
number  expressions  (16)  and  (17)  gives  the  following: 


/l  + a - 1 


-C  /I  + a 


/l  - a 


l - /n 


♦ (a  < l ! 


-C  vT  - a , (a  < 1 ) . 


The  mean  energy  for  each  wave  may  be  expressed  as  follows: 

E=tp9A2  (33) 

n 6 n’ 


where  the  amplitudes  of  the  waves  A^ , A^  and  A^  are  taken  from  (26),  (27), 
(30)  and  (18).  Substituting  the  group  velocities  (32)  and  the  energies  (33) 
into  the  work  rate  expression  (31)  gives  the  following: 


f («), 
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where  f(a)  * f^(a)  + f2(a)  and 


*1<«>  "f  0 + f> 


f2(a)  = 


#1  3 v 

'2  ~ 8 

/l  - a~ 


* a < 1 
, a > 1 


The  normalized  work  rate  function  f(a)  is  shown  in  Figure  5.  One  sees  the 
familiar  resonance  at  a = 1.  The  resonance  occurs  because  the  work  performed 
by  the  pressure  distribution  on  the  fluid  creates  energy  that  cannot  propagate 
away  from  the  disturbance.  This  follows  from  equation  (32)  where  one  sees 
that  the  relative  group  velocities  of  the  third  and  fourth  waves  are  zero 
for  a = 1.  The  expressions  for  the  group  velocities  (32)  and  the  mean  work 
rate  (34)  are  consistent  with  Wu's  [5]  results. 

THE  SURFACE  ELEVATION  FOR  THE  UNIFORM  PRESSURE  DISTRIBUTION 
The  surface  elevation  Z(x,t)  due  to  a uniform  pressure  distribution  was 
computed  from  c(x,t)  using  the  superposition  integral  (8).  The  surface 
elevation  caused  by  the  delta  function  distribution  c(x,t)  was  obtained  from 
(18),  (24),  (26),  (22)  and  (27)  - (30).  One  may  write  the  uniform  pressure 
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The  results  are  expressed  in  nondimensional  form  as  follows: 
z'  = — = (a  + ib)  e;lot  (36) 

P0/Dq 

where  z'  is  the  normalized  wave  elevation,  a is  the  component  in-phase  with 
the  pressure  and  b is  the  out-of-phase  component.  Expressions  for  the  in  and 
out"of-phase  components  for  a < 1 are  given  in  the  Appendix.  The  two 
components  of  the  wave  elevation  are  functions  of  the  normalized  frequency 
a and  the  Froude  number 

F = c/(g  l)y2m 

Numerical  results  for  F=0.7  are  shown  in  Figures  6,  7 and  8 for  three 
frequencies:  a=0,  0.5  and  0.95.  In  each  figure  the  normalized  in-phase  (a) 
and  out-of-phase  components  (b)  are  plotted  as  a function  of  the  nondimensional 
distance  x'  = x/l.  Figure  5 shows  the  in-phase  component  of  the  wave 
elevation  for  the  zero  frequency  case.  (The  out-of-phase  component  is  zero). 
The  near  field  disturbance  resembles  the  wake  produced  by  a planing  surface. 
The  standing  wave  in  the  far  field  downstream  is  evident.  Figures  6 and  7 
show  the  in-phase  and  out-of-phase  components  of  the  wave  elevation  for 
a=0 . 5 and  0.95,  respectively.  One  can  see  interference  between  the  various 
waves  in  the  far  field  downstream  for  both  frequencies.  In  addition,  a long 
wave  appears  upstream  for  a=0.95,  but  is  not  apparent  for  a=0.5. 

To  clarify  the  behavior  of  the  waves  in  the  far  field,  Table  1 has  been 
prepared.  The  ratio  of  wavelength  to  pressure  distribution  length  (x/£)  has 
been  calculated  for  each  wave  at  all  three  frequencies  using  the  wavenumber 
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expressions  (16)  and  (17).  One  sees  from  the  table  that  the  second  and 
fourth  wave  have  the  same  wavelength  and  form  a standing  wave  for  a=0, 
because  one  wave  travels  to  the  right  and  the  other  to  the  left.  The  first 
and  third  waves  have  infinite  wavelength,  but  their  amplitudes  are  zero. 

For  nonzero  frequencies,  the  length  of  the  second  wave  is  shorter  than 
the  fourth.  This  causes  the  interference  pattern  in  the  downstream  wake. 
The  third  wave  occurs  upstream.  Its  amplitude  is  too  small  to  appear  in 
Figure  6,  but  it  is  evident  in  Figure  7.  The  wave  is  so  long  at  a=0.95 
that  only  about  half  a cycle  of  the  wave  appears  in  the  Figure. 

TABLE  1 

WAVELENGTHS  IN  FAR  FIELD  PRODUCED  BY  UNIFORM  PRESSURE  DISTRIBUTION 


a 

A z/l 

V* 

a4/£ 

0 

oo 

3.08 

CO 

3.08 

0.5 

244 

2.49 

1 .44 

4.23 

0.95 

78.3 

2.14 

20.4 

4.93 

SWWARY  OF  RESULTS 

The  steady  state  surface  elevation  is  expressed  in  (18)  as  the  sum  of 
four  integral  terms,  each  having  a simple  pole  singularity  in  the  complex 
wavenumber  plane.  The  location  of  the  poles  is  shown  in  Figure  2.  When 
the  Stokes  viscosity  is  suppressed  all  the  poles  lie  on  the  positive  real 
axis  for  low  frequencies  (a<l).  At  higher  frequencies  (a>l)  only  the  first 
two  lie  on  the  real  axis.  The  integration  paths  are  indented  as  shown  in 
Figure  3 so  that  the  poles  lie  on  the  same  side  of  the  path  as  they  did 
with  Stokes  viscosity.  Each  integral  term  is  then  interpreted  as  a path 
integral  in  the  complex  wavenumber  space  as  indicated  in  equation  (21). 
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To  reduce  the  number  of  integrations,  three  symmetry  relations  that  follow 
from  the  integral  forms  and  the  integration  paths  are  introduced. 

The  first  integral  form  is  evaluated  by  selecting  appropriate  closed 
contours  for  the  positive  and  negative  x-  axis,  as  shown  in  Figure  4,  and  then 
applying  Cauchy's  residue  theorem.  The  integral  is  expressed  in  discontinuous 
form  in  equations  (24)  and  (26).  The  integral  consists  partially  of 
auxiliary  exponential  integrals  which  are  known  tabulated  functions.  In 
addition  a residue  term  appears  on  the  negative  x-axis  or  the  downstream  side. 
This  term  is  an  exponential  and  it  represents  the  familiar  undamped  surface 
wave.  The  other  three  integrals  are  written  from  the  first  result  using 
the  symmetry  relations. 

The  surface  elevation  is  seen  to  consist  of  exponential  integral  terms 
contributing  to  the  near  field  and  exponentials  appearing  on  only  one  side 
of  the  disturbance.  The  exponential  terms  are  the  free  waves  and  are 
associated  with  the  residues  of  each  pole.  For  low  frequencies  (a  < 1 ) 
three  waves  appear  downstream  and  one  upstream,  as  has  been  noted  by  previous 
investi gators . [3,4,5] 

At  the  higher  frequencies  two  of  the  poles  leave  the  real  axis.  The 
exponential  integral  terms  in  the  solution  remain  the  same,  except  that  their 
arguments  become  complex  instead  of  real.  The  residue  terms  persist  until 
a=2 , at  which  point  the  poles  leave  the  integration  contour.  For  1 < a < 2 
the  residue  terms  decay  exponentially  with  distance  and,  therefore,  do  not 
contribute  to  the  far  field  solution. 

The  mean  rate  at  which  energy  propagates  away  from  the  disturbance 
was  computed.  The  result,  which  agrees  with  Wu  [5],  is  given  in  equation  (34), 
and  is  shown  in  Figure  5. 
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The  surface  elevation  was  calculated  for  a uniform  pressure  distribution. 

This  was  done  to  eliminate  the  logarithmic  singularities  that  appear  in  the 
solution  for  the  delta  function  distribution.  Expressions  for  the  normalized 
in-phase  and  out-of-phase  components  of  the  surface  elevation  are  given  in  the 
appendix  for  the  low  frequency  case.  Numerical  results  are  presented  in 
Figures  6,  7,  and  8 for  a Froude  number  of  0.7  and  u=0,  0.5,  and  0.95, 
respectively.  The  zero  frequency  result  (a=0)  shows  that  the  water  surface  deforms 
like  that  of  a planing  surface  in  the  near  field,  while  standing  waves  are 
apparent  in  the  far  field.  For  the  nonzero  frequencies  (^=0.5  and  0.95)  inter- 
ference occurs  in  the  waves  downstream.  In  addition,  a long  wave  occurs  upstream 
for  t=0.95. 

The  lengths  of  the  various  waves  were  calculated  for  the  three  frequencies. 

The  results  are  shown  in  the  table.  The  zero  frequency  standing  wave  is  shown 
to  consist  of  two  downstream  waves,  each  having  the  same  wavelength  and 
traveling  in  opposite  directions.  The  interference  pattern  in  the  downstream 
wave  pattern  for  nonzero  frequency  is  caused  primarily  by  the  same  two  waves, 
whose  wavelengths  now  differ.  A long  wave  appears  upstream  as  one  approaches 
the  critical  frequency. 
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APPENDIX 


EXPRESSIONS  FOR  IN  PHASE  AND  OUT  OF  PHASE  COMPONENTS  OF 
SURFACE  ELEVATION  FOR  UNIFORM  PRESSURE  DISTRIBUTION 


The  normalized  in  phase  component  of  the  surface  elevation  is  given 
as  a(F,a)  and  the  out  of  phase  component  as  b(F,a).  Each  component 
consists  of  four  parts  corresponding  to  the  four  poles  , K^,  and  K^. 
One  may  write  for  a and  b: 


a = 


a . 
1 


and 


4 


V“ 


b = 


b.  . 


i'  = i 1 

As  a result  of  the  discontinuous  representation  of  the  surface  elvation 
C(x,t),  the  superposition  integration  (8)  must  be  performed  in  three 
regions:  1)  upstream  (x’>  1/2),  2)  Tinder  the  pressure  distribution  (Ixl  < 1/2) 
and  3)  downstream  (x ’<  - 1/2). 

The  arguments  for  the  exponential  integral  and  exponential  terms  are 
given  in  nondimensional  form.  First,  the  wavenumbers  are  normalized  as  follow??: 


y . (a)  = K./K  , i = 1,2, 3,4. 
l l o 

1/2 

The  speed  is  expressed  in  terms  of  the  Froude  number  F = c/(g£)  , and 

the  normalized  longitudinal  distance  as  x ' = 2x/£. 
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b~  = < g(Y  1^)  - g(Y«L  ) ♦ lnW^L  ) - ln(y  L ) 

2ir/  i+a  ^ 1 1 1 2 1 1 J 


■*3  = ^ ^ jg(Y3Ll>  ‘ g(Y3L2)  + ln(Y3Ll)  ' ln(Y3L2)  + 


- 4 it  cos  (v3X)  sin  (y ^G) 


,9(Y4L1)  - 9(Y4L2)  + InfY^)  - ^ (>41-2) 


|x'|<  1/2  (Under  Pressure  Distribution),  a<l 


a.  = < tt  + f(Y  L ) + f (Y  L ) -2Y  cos  (Y  L ) 

1 2ir^i+«  13  11  13 


Zl j*  + f(72L3)  + 

2tr/  1+a  : 


ffY2L3)  + f(Y2Ll)  ~2tt  cos(V2L3)> 


- * <"  + f(Y-L..)  + fa,L.)  ♦ *2tt  cos  (y  L ) 

3 2^/i^r  ] 3 3 3 1 3 1 


j = ^ + f(y  Ln)  - ffy^L,  ) - 2tt  cos  (y  ,L.J 

4 ] 4 3 4 1 4 3 


b = — — M-iV  " 9^IL3>  + In  (y  L > - In  (y^) 

+ 2tt  sin  ^Y^L3)t 


o2  = Sg(Y2Ll)  " 9(Y2L3)  + ln  (Y2L1)  " ln  (Y2L2) 

2-n  /l+a  ^ 

+ 2tt  9in  CY  2 L3 ) . 
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pn^PH — - 


r 


-1 


2tt/  1-a 


^g(Y3L3)  - 9(Y3L3)  + ln(Y3L3)  - lnty^) 
+ 2n  stn(Y^L^)> 


b = — r-  <(g(Y4L3)  - 9(7^)  ♦ ln<Y4L3>  - 1"(^4) 

2tt/  1-a  C > 


- 2-n  sin(Y.L  )> 

4 4 


3.  x #<  - 1/2  (Downstream),  a<l 


-1 


3i  = — ^ f (yxl4)  - f(y1L3^  + 4ir  sin  (y^x)  sin  (^iGv 


= 


2ttv  1+a 
1 


2 2ir/  1+a" 


-1 


3 2*/^r 


a . = 


4 2tt/  l-a_ 


b,  = 


-1 


1 2tt/  l+"a~ 


b2  = 


2tt/  1+a 


b3=  " 


2tt/  1-  a 
4 2tt/  1-  a 


<^f(Y2L4)  - f(Y2L3}  + 4lT  sin  (Y2X)  Sin  (Y2G)^ 


<f(Y3V  - f(V3^ 


<,f(Y4L4)  - f(Y4Lj)  + 4 it  sin  (Y4X)  sin  (y4g| 


|g(Y1L4)  - g(Y1L3)  + lnfY^)  - In  (y^j)  + 
4 tt  COS  (Y  X)  sin  ( Y x G ) i 


g(Y2L4>  - g(Y2L3)  + In  (y2l4)  - In  (y2l3)  + 


4 it  cos  (Y2X)  sin  (Y2G)j 


9(Y3L4)  " 9<Y3LS)  + lnt73L4>  ' ln(Y3L3) 


9<Y4V  - 9<V3}  + ln(Y3V  " ln<Y4L3)  + 


4tt  cos  (Y4X)  sin  (“Y 4G^J 
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a.  Migration  of  K,  and  K 
1 2 


Figure  2 - Migration  of  Poles  in  Complex  Wavenumber  Plane  as 
Frequency  Parameter  a Increases 


r 


K, 


r, 


a.  Path  for  I, 


K, 


b.  Path  for 


K3 

o- — & — -r. 


c.  Path  for  1^,  a < 1 


d.  Path  for  a < l 


K- 


* 


r3*r4 


K. 


e.  Path  for  I and  I for  a > 1 

-5  4 


Figure  3 - Integration  Paths  in  Complex  Wavenumber  Plane  for 
Dissipationless  Liquid 
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Figure  5 - Normalized  Mean  Work  Rate  Function  f(a)  as  a Function 
of  the  Frequency  Parameter  a 
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a.  .In  phase  component 


b.  Out  of  phase  component 


Figure  7 - Normalized  Wave  Elevation  for  F=0.7  and  a=0.5 
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"V 7: 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


(1)  DTNSRDC  REPORTS.  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF  ] 

PERMANENT  TECHNICAL  VALUE.  DESIGNATED  BY  A SERIAL  REPORT  NUMBER  j 

(2)  DEPARTMENTAL  REPORTS.  A SEMlFORMAL  SERIES.  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE.  OR  OF  LIMITED  INTEREST  OP 
SIGNIFICANCE.  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION 

(3)  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS.  NUMBERED  AS  TM  SERIES 
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